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ABSTRACT 

We introduce a nonextensive entropic measure S x that grows like N x , where N is the size of the system 
under consideration. This kind of nonextensivity arises in a natural way in some iV-body systems endowed 
with long-range interactions described by r~ a interparticle potentials. The power law (weakly nonextensive) 
behavior exhibited by S x is intermediate between (1) the linear (extensive) regime characterizing the standard 
Boltzmann-Gibbs entropy and the (2) the exponential law (strongly nonextensive) behavior associated with 
the Tsallis generalized g-entropies. The functional S x is parametrized by the real number ^ £ [1, 2] in such a 
way that the standard logarithmic entropy is recovered when \ = 1 . We study the mathematical properties 
of the new entropy, showing that the basic requirements for a well behaved entropy functional are verified, 
i.e., S x possesses the usual properties of positivity, equiprobability, concavity and irreversibility and verifies 
Khinchin axioms except the one related to additivity since S x is nonextensive. For 1 < % < 2, the entropy 
S x becomes superadditive in the thermodynamic limit. The present formalism is illustrated by a numerical 
study of the thermodynamic scaling laws of a ferromagnetic Ising model with long-range interactions. 

Keywords: entropy, nonextensivity, long-range interactions. 
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1. INTRODUCTION 



There is nowadays an intense research activity on the mathematical properties and physical applications 
of new versions of the Maximum Entropy Principle based on generalized or alternative entropic functionals 
|j] |l8|. This line of inquiry has been greatly stimulated by the work of Tsallis, who showed that it is 
possible to build up a mathematically consistent and physically meaningful generalization of the standard 
Boltzmann-Gibbs-Jaynes thermostatistical formalism on the basis of a nonextensive entropic measure 0. 
The main motivation behind Tsallis proposal is that there are important physical scenarios, such as self- 
gravitating systems fl9| , [20|| , electron- plasma two dimensional turbulence among many others, which 
are characterized by a nonextensive behavior: due to the long range of the relevant interactions some of 
the thermodynamic variables usually regarded as additive, such as the internal energy, lose their extensive 
character. This suggests that a nonextensive (non-additive) entropy functional might be appropriate for 
their thermostatistical description. The Jaynes MaxEnt approach to Statistical Mechanics J2^,|23| suggests 
in a natural way the possibility of incorporating alternative entropy functionals to the variational principle. 
The new entropy functional introduced by Tsallis has the form 

where (pi, i = 1, . . . ,w) are the microstate probabilities describing the system under consideration and the 
entropic index q is any real number. The standard Boltzmann-Gibbs entropy S = — yi_i Pi In Pi is recovered 
in the limit q — > 1. The measure S q is nonextensive: The entropy of a composite system A © B constituted 
by two subsystems A and B, independent in the sense that p\f® B%> — p\ Pj , verifies the Tsallis' g-additive 
relation 

S q (A ®B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B). (2) 

We see from the above equation that the Tsallis' parameter q can be regarded as a measure of the degree of 
nonextensivity. Many relevant mathematical properties of the standard thermostatistics are verified by the 
Tsallis' generalized formalism or can be appropriately generalized p]JIq|. Self-gravitating systems constituted 
the first physical problem discussed within Tsallis' nonextensive thermostatistics and Tsallis' theory has 
recently been applied to other physical problems j24|-|2{| . One of the main consequences of the intensive effort 
devoted in recent years to the study of Tsallis theory is that there is now a growing consensus that there 
are many problems in statistical physics, biology, economics, and other areas, where a generalization of the 
standard approach based on Boltzmann-Gibbs-Jaynes extensive entropy might be useful. A comprehensive 
bibliography on the current research literature on Tsallis' theory and the statistical physics of nonextensive 
systems can be found in ||(]|3l| . 

Inspired on Tsallis' pioneering proposal, various nonextensive entropic measures endowed with interesting 
properties have been recently advanced ]q-p^]. Moreover, it has been proved that some physically relevant 
mathematical properties are shared by large families of entropic measures . The aim of the present work 
is to explore the possibility of developing a thermostatistical formalism based on a nonadditivc entropic 
functional characterized by a degree of nonextensivity weaker than the one exhibited by Tsallis measure. 
As we shall see, Tsallis entropy varies exponentially with the size of the system under consideration. The 
new measure introduced here only scales as a power of the size of the system. That is to say, our proposal 
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is tantamount of considering a nonextensive regime intermediate between (1) the standard extensive one 
associated with Boltzmann-Gibbs entropy, and (2) the exponentially nonextensive one described by Tsallis 
formalism. 

There is an important physical motivation for introducing an entropy endowed with power law nonextensiv- 
ity Systems with long range interactions constitute the physical scenarios where the need of a generalization 
of the standard thermostatistical formalism can be more clearly appreciated. Consider a system of TV par- 
ticles in a d-dimensional (one particle) configuration space. If the dependence of the interparticle potential 
with the interparticle distance r is given by ?*~ Q , it can be shown that the system's energy levels scale as |32| 

TV TV — TV £!- (3) 

1 — a/a 

Hence, for large TV, the internal energy scales as a power of the size of the system. That is, 

E ~ N 2 - a/d (4) 

In the case of extensive systems governed by short-range interactions, the internal energy E and the the 
entropy S scale in the same way: they both grow linearly with TV. On the other hand, the temperature 
T is an intensive variable and does not change with TV. How can these scaling laws be generalized to the 
non-extensive setting? A possible path towards the alluded generalization starts with the Helmholtz free 
energy, 

F = E - TS. (5) 

From the above expression it seems reasonable to expect both E and TS to scale in the same fashion p2| . 
The only way to fulfill this expectation, if the standard extensive entropy is used, is to require that the 
temperature scales as [ f32"|]33"| 

T ~TV ~ N 1 ~ a / d j (6) 

losing its intensive character. 

It would be very appealing to have, within the nonextensive scenario, an entropy S x endowed with the 
same scaling law as the one exhibited by the energy. If, as it occurs for extensive systems, the entropy and 
the energy behave in the same way, the temperature would still be an intensive quantity. What we need in 
order to have a thermostatistical formalism with these nice properties is an entropy functional with power 
law nonextensivity, scaling as 

s x ~ m ~ N 2 - a ' d . (7) 

From the above discussion it is clear that we are going to assume that the exponent x appearing in the 
cntropic scaling law is related with d and a by 

X = 2 - a/d, (8) 

so that the physically motivated range of vales for \ i s 

1 < X < 2. (9) 
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The purpose of this work is to study the basic properties of a possible candidate for a weakly-nonextensive 
thcrmostatistics, and to consider its application to a magnetic system with long-range interactions. The 
paper is organized as follows. In Section II the exponential behavior of Tsallis entropy is analyzed and 
a power-law weakly nonextensive entropy is introduced. The basic mathematical properties of the new 
measure are studied in Section III. Two state systems are considered in Section IV. In Section V the weakly 
nonextensive entropy is applied to a ferromagnetic Ising model with long-range interactions. Finally, some 
conclusions are drawn in Section VI. 



2. EXPONENTIAL VS. POWER LAW NONEXTENSIVITY 

Tsallis Entropy and Exponential Nonextensivity 

The nonextensive behavior of Tsallis functional, encapsulated in equation (|^), is the most important 
single property distinguishing Tsallis measure from the standard additive logarithmic entropy. An important 
consequence of relation (||) that has not been fully appreciated in the literature is that Tsallis entropy may 
varies exponentially with the size of the system. In order to clarify the above assertion let us consider a 
system consisting of N independent two-state subsystems. For simplicity we also assume that the system is 
described by an equiprobability distribution. That is, each of the w = 2 N possible microstate of the system 
has probability 1/w. The associated Tsallis measure adopts then the value S% quiproh - = 1- ™ * , so that, for 
q < 1 and large N, 

S q — 2 (1 -^ N . (10) 

i - q 

If q > 1 the g-entropy tends to the constant limit value — ^j- as N — > oo. We are going to consider only the 
q < 1 regime. In that case the g-entropy exhibits an exponential behavior as a function of the size N of the 
system. In general, the entropy S q (N) of a composite system A = ®^ =1 A^ consisting on N identical and 
independent subsystems {A^\ j = 1, . . . N}, is given by 

S g (N) = -1_ {1 - [1 + (1 - q)S q (l)] N } , (11) 

where S q (l) stands for the individual entropy of each constituent subsystem A^>. Notice that in order 
to obtain the above expression it is not necessary to assume that each subsystem A^) is described by an 
equiprobability distribution. 



Power Law Nonextensive Entropy 

We believe that it is worth exploring the possibility of a thermostatistical formalism based on a nonex- 
tensive entropy S x growing as a power N x of the size of the system, instead of growing exponentially. In 
accord with the physical arguments discussed in the Introduction (see equations (|8||j|)), we shall restrict our 
considerations to values of \ £ [1, 2] . Let us assume that the functional S x associated with a given discrete 
probability distribution {p,;, i = 1, . . . , w} is given by 

w 

S x = Y,<t ) x(Pi), (12) 

i=l 
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4> x (pi) being an appropriate function of the individual microstate probabilities pi. In order to find a suitable 
expression for the function <fi x (z) it is enough to consider again the equiprobability distribution associated 
with a collection of N identical independent two state subsystems. In that case we have w = 2 N and 



S x = w <t>x( l l w ) 

= 2 N Mz- N ) 



(13) 



which implies 



,(2-*) ~ 2~ N N X . 



(14) 



The simplest choice for a function (f> x (z) complying with the above relation is 

<^x( z ) = z(-\nz) x . 



(15) 



Unfortunately, this function is not adequate for our purposes. Its second derivative d 2 (f> x /dz 2 lacks a definite 
sign within the relevant range of values of z and x, leading thus to an entropy functional without a definite 
concavity. However, since we are only interested in the large N asymptotic regime, any function <Px( z ) 
behaving like (15) in the limit z — > would do. As we shall see, the function 



leads to the measure 



,(z) = - z(-\nz) x + z(-lnz) 1/x 



1 w 

S * = 2 ^ [pi{-^Pi) x +Pi{-^Pi) 1/x 



(16) 



(17) 



which complies with all the basic properties of a well behaved entropy. It is clear that in the case x = 1 the 

w 

standard entropy Si = — V^Pi lnpj is recovered. 

i=l 

The optimization of S x under the constraints imposed by normalization, 



5^ Pi = i, 



(18) 



and the mean values 



(AU) = Y.P* A i = A x T) (r = l,...,Ji) 



(19) 



of a given set {^4*- r ' } of observable leads to the variational problem 



S \ S x- X>M (r) > -aEl" 



0. 



r=l 



whose solution is of the form 



(20) 



Pi = F\a + Y j (3 r A 



(r) 



t = l, 



(21) 



Here a and {/3 r } are appropriate Lagrange multipliers and F(z) is the inverse function of 
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(— In z) x + (— In z) x — %(— In z) x 1 (— In z) x ' 

X 



(22) 



which verifies the relation 



FW x {z)] = 4>' x [F(z)\ = z. 
The function F{z) is well defined because the second derivative of s x (z) is negative, 



dz 2 



2z 



X (x-l)(-lnz) x - 2 +1(1-1] (-lnz)t 
X \X / 



— %(— lnz) x 1 H — (— InzW 

X 

<0, 



(23) 



(24) 



for < z < 1, and 1 < x < 2. Actually, s^(z) has a definite (negative) sign within the larger interval 
[1/XcjXcL where Xc ~ 2.1762. This can be appreciated in Fig. [I], where the real roots of s' x (z) = are 
depicted as a function of x- We see that there are no real roots \ G [1/xc Xc]- 



z 




0.125 0.25 0.5 1 2 4 8 



X 

FIG. 1. The real roots of s"(z) = as a function of \. There are no real roots within the interval [l/Xc,Xc], 
where X c ~ 2.1762. 

It is important to realize that the lack of a simple analytical expression for the function F(z) does not 
constitute a serious conceptual problem, nor does it pose any practical difficulty for the numerical treatment 
of problems involving the MaxEnt distributions (^l|). The form of the function F(z) is shown in Fig. ||. 



G 




FIG. 2. The function F(z) denned by the equations (g2jj23|) for x = 1-0, 1.25, 1.5, 1.75, 2.0. 

The canonical ensemble probability distribution associated with the entropy S x , 

Pi = F(a + (3a), (25) 

is obtained when the S x is extremalized under the constraints of normalization and the mean value of the 
energy, 

(U) = £>e;, (26) 



where q stands for the energy of the microstate i. The Lagrange multiplier j3 appearing in (|25|) is the one 
associated with the energy constraint and corresponds, within the present thermostatistical formalism, to 
the inverse temperature. That is, taking Boltzmann constant k = 1, we have f3 = l/T. 

3. MAIN PROPERTIES OF THE NEW ENTROPY 
3.1 Khinchin axioms 

Khinchin proposed a set of four axioms p4[ , which are usually regarded as reasonable requirements for 
a well behaved information measure. Our entropy measure S x verifies the first three of them: 

(i) S x = S x (pi, . . . ,p w ), 

i.e., the entropy is a function of the probabilities pi only. 

(ii) S x ( Pl ,..., Pw ) < S X (A,...,A) = S^-(w), 

i.e., S x adopts its extreme at equiprobability (this property will be proved in subsection 3.2). 

(iii) S x (p!,...,p w ) = S x (pi,...,p w ,0), 

this property, known as expansibility, is clearly verified since 4> x (0) = 0. 



7 



(iv) The fourth Khinchin axiom concerns the behavior of the entropy of a composite system in connection 
to the entropies of the subsystems. We will comment on this axiom later. 



3.2 General mathematical properties 

Let us now consider the properties related to positivity, certainty, concavity, equiprobability, additivity 
and irreversibility, which are some of the most important features characterizing an information or entropic 
measure ifflliJIlfU 



3.2.1 Positivity 

It is plain from equation ( |l6| ) that 4> x (0) — 4>xW = an d also that <f>x(p) > for p G [0, 1]. This implies 
the positivity condition 

S x > 0. (27) 



3.2.2 Certainty 

The equality symbol in Eq. (I^ ) holds only at certainty, i.e., 

S X (1,0,...,0) = 0. (28) 

Indeed, S x vanishes if and only if we have certainty. 

3.2.3 Concavity 

Considering (pi, . . . ,p w ) as independent variables, the second partial derivatives of S x are 

pa a f)2o 

x - x 5 jk < for < Pj < 1. (29) 



dpjdpk d 2 Pj 

Expression (|29| ) guarantees definite concavity over probability space (see, for instance, fj3f). 

3.2.4 Equiprobability 

The probability distribution that extremizes S x under the normalization constraint is 

Pi = F(a) = 1/w (30) 

Therefore, since S x has negative concavity, it is maximal at equiprobability. A well behaved entropy should 
also be, at equiprobability, a monotonically increasing function of the number of states w. For large w we 
have S x (w) ~ (lnu>) x . Therefore, for large values of w, S x qutpr ' (w) is an increasing function of w. 
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3.2.5 Nonextensivity 



The nonextensive behavior of the entropy S\ is determined by the relation of the entropy of a composite 
system with the individual entropies of its constituent subsystems. Let us consider systems A and B with 
associated probabilities {a^, i = \, . . . ,wa\ and {bj, j — 1, . . . ,wb}, respectively. If systems A and B 
are independent, i.e., the composite system A © B has associated probabilities {a^; i — 1, . . . ,wa', j = 
1, . . . , wb}, then the entropy S X (A © B) of the composite system minus the sum of the entropies of its 
subsystems, 

AS X (A, B) = S X (A (SB)- S X (A) - S X (B) (31) 

is the quantity characterizing the nonextensive features of the measure S x . When AS X (A, B) > 
0(AS x (A,B) < 0) we have superaditivity (subadditivity) . From the examination of particular examples 
we conclude that AS X (A, B) does not have always the same sign. However, the region of probability space 
where AS X (A, B) is positive is much larger than the region where that quantity is negative. Furthermore, 
if we consider N identical subsystems (instead of just two of them), the region of probability space cor- 
responding to subadditive behavior tends to vanish as N grows, consequently, the entropy S x becomes 
superextensive in the thermodynamic limit. Particular examples illustrating these features of the measure 
S x are provided in Section IV. Moreover, in Section V we are going to consider the nonextensivity of S x in 
connection with the thermodynamic properties of an Ising model endowed with long-range interactions. 



3.2.6 Irreversibility 

One of the most important roles played by entropic functionals within theoretical physics is to characterize 
the "arrow of time" . When they verify an ii-theorem, they provide a quantitative measure of macroscopic 
irreversibility. We will now show, for some simple systems, that the present measure S x satisfies an H- 
theorem, i.e., its time derivative has a definite sign. 

Let us calculate the time derivative of S x 

f = i <»> 

i— 1 

for a system whose probabilities pi evolve according to the master equation 

-§ - YpXPi-PiiPi], (33) 

where Pij is the transition probability per unit time between microscopic configurations i and j. Assuming 
a system with a uniform equilibrium distribution and detailed balance, i.e., Pij = Pji, we obtain from ( |32| ) 

dS x 



= \ E E P ^ (W - ft) ( s xfe) - s 'xW ■ (34) 

4=1 j = l ^ ' 

Since s''(j>i) < 0, the quantities (pi — Pj) and (s' x (pj) — s' (pi)) have the same sign. Then we obtain 
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The equality holds for equiprobability, i.e., at equilibrium, while in any other cases the entropy S x increases 
with time. Therefore, S x exhibits irreversibility. 



3.2.7 Jaynes thermodynamic relations 



It is noteworthy that, within the present MaxEnt formalism, the usual thcrmodynamical relations involving 
the entropy, the relevant mean values, and the associated Lagrange multipliers, i.e., 

^% - * < 36 > 

are verified. Hence, our formalism exhibits the usual thermodynamical Legendre transform structure. Actu- 
ally, this property is verified by a wide family of entropy functionals [^Q. A particular important example 
of (|3^) is furnished by the canonical ensemble thermodynamic relation 

- - = P = 7F- ( 37 ) 



d(U) ^ T 



4. TWO-STATE SYSTEMS 



In order to illustrate some of the above properties, we consider a two-state system (with associated 
probabilities {p, 1 — p}). In this case, S x only depends on the variable p. In fact, from its definition, we have 



S x(p) = o P(- ln P) X +p(-mp)* 

+ (1 - p)(- ln(l - p))* + (1 - p){- ln(l - p))- 



(38) 




FIG. 3. Entropy S x for a two-state system (w = 2) as a function of p for different values of \ indicated on the 
figure. 
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The shape of S x (p) for different values of \ is shown in Fig. |3|, which exhibits the positivity and concavity 
of S x . In fact, from expression (|38|), the first derivative of S x (p) vanishes at p = 1/2 and d 2 S x /dp 2 < Vp G 
[0,1]. Since the second derivative is always negative, S x (p) is maximal at equiprobability. Moreover, as 
shown in the general case, taking into account the concavity of S x and that S x vanishes at the certainty, 
then S x is positive for all p. 

The non-additivity of S x is illustrated in Fig. |J for two independent two-state systems A (with probabilities 
p and 1 — p) and B (with probabilities q and 1 — q), through the plot of the relative difference 

(AS x ) reL = [S X (A ® B) — (S X (A) + S X (B))]/S X (A © B) (39) 

as a function of p and q. For most values of p and q the nonextensive measure S x behaves in a superad- 
ditive fashion. Only for values lying in a small region near the edges of the (p, q)— square does S x become 
subextensive. 



0.01 - 




FIG. 4. The subadditivity and superaditivity of S x . The quantity (AS x ) r ei. as a function of p and q, for a 
composite system constituted by the independent two-state subsystems A (with associated probabilities {q, 1 — q}), 
and B (with associated probabilities {p, 1 — p}). 



The nonextensive behavior of S x in the thermodynamic limit can be illustrate by recourse to a system 
constituted by N two-state subsystem. Let us first assume that each one of the subsystems is described by 
the same probabilities p and (1 — p). Consequently, the entropy associated with the composite system is 
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1 N ( N\ r 



It can be shown after some algebra that, for < p < 1, 

[-ln(l -p)] x N x < S X (N) < [-lnp} x N x + [- lnp] 1/x iV 1/x . 
Hence, for any given value of p € (0, 1/2) there exist an M such that 

N > M ==> 5 x (iV) > NS X (1), 



(40) 



(41) 



(42) 



which means that Sx becomes superadditive for large enough values of N. This is shown in Fig. ||a , where 
the quantity 



(AS X U. = [S X (N) - NS X (1)]/S X (N) 



(43) 



is depicted for different values of N. In a similar way, Fig. ||b shows the behavior of (AS x ) re i. as a function 
of p for a composite system consisting of TV — 1 two-state subsystems each with probabilities {5,5} along 
with one extra subsystem with probabilities p and (1 — p). 



in 



0.06 



0.03 



-0.03 



N=30 


(a) 


N=3 











0.25 



0.75 



on 



0.25 



0.06 - 


N=30 


(b) 




N=3 




0.03 






- 


N=2 








0.03 







0.5 
p 



0.75 



12 



FIG. 5. The superaditivity of S x in the thermodynamic limit, (a) The quantity (AS x ) re i. associated with a 
system constituted by -/V identical two-state subsystems with probabilities p and (1 — p), as a function of p and for 
different values of N. (b) (AS x ) re i. as a function of p for a system constituted by N— 1 identical two-state subsystems 
with probabilities {i, i } along with one extra two-state subsystem with probabilities p and (1 — p). 



5. FERROMAGNETIC ISING MODEL WITH LONG-RANGE INTERACTIONS 



In this section we apply the canonical thermostatistics associated with the entropy S x to an ./V-body system 
described by an r~ a intcrparticlc potential. The energy of this kind of systems scales as a power of the 
number of particles ]32| . The consideration of systems exhibiting this power-law nonextensivity constituted 
our main motivation for introducing the measure S x . As argued in Section I, an entropy endowed with 
the same non-extensive behavior as the one associated with the energy may lead to a thermostatistical 
description preserving the intensive character of the temperature. In order to illustrate these ideas we are 
going to study a long-range Ising model described by the Hamiltonian 

H = J ^ 1 -4^ > (Si = ±l,Vi), (44) 

i,j=l 

where J is an appropriate coupling constant, and the sum runs over all the pairs of sites on a d-dimensional 
lattice with periodic boundary conditions and stands for the distance between the sites i and j. It is clear 
that the range of the interaction is determined by the value of the exponent a. In particular, the standard 
(short-range) first-neighbors interaction is recovered in the limit a — > oo, while the mean field approximation 
is obtained when a = (replacing J by J/N). These extreme cases illustrate, respectively, the two possible 
thermodynamic behaviors that, according to the value of a, are admitted by the system (Q). On the one 
hand we have the extensive regime corresponding to a > d. On the other hand we have the non-extensive 
regime associated with a < d. In order to clarify this let us estimate the internal energy per particle at zero 
temperature. We have Q 



E 
N 



drr d - l r- a . (45) 



It is easy to see that the above integral converges if a > d and diverges if < a < d. The following is a 
standard notation used in previous works: 



1 f N ., , N^^-a/d 

N = l + -I drr d - 1 r- a = - (46) 

d Ji 1 — a/d 

Here we consider the non-extensive regime for the one-dimensional case d = 1 and a = 0.8. We have 
performed numerical simulations by using a novel approach recently introduced to study systems with 
long range interactions governed by generalized entropies as the one considered here. The method relies 
upon the calculation of the number of states, fl(ek), with a given energy efe. 

Note that, the number of possible configurations and associated probabilities pi is in general very large, 
W = 2 N for Ising models for example. However the number of permitted energies or energy levels K is not 
so large, because there is a large number of states tt(ek) sharing the same energy ek- We can rewrite the 
sums in Eqs. (|l7 26) taking into account the K energy levels weighted by fi(efc): 
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1 K 

S x = -^n( £fc )p(e fc )[(-ln(p(e fc )) x + (-ln(p(e fc )) 1/x ], (47) 
fe=i 
if 

<0> =5^n(e fc )p(efc)0(e fc ). (48) 
fc=i 

Hence, the knowledge of the number of states fi(efe) allows, by the use of these expressions, the calculation 
of the entropy and any averages of interest. 

To calculate the probabilities p(ek) we use the definition p(e&) = F(a + /3efc), where the function F is 
compute numerically as the inverse function of s'(z) (See Fig. ||). 

To compute the histogram fi(efc), the "Histogram by Overlapping Windows (HOW)" method [p5Pj36t 
is used. A naive way of computing the histogram consists in generating different system configurations 
randomly and counting how many times a configuration with energy ej, appears. However, since the f2(efc) 
values span too many orders of magnitude it is not possible to find in this way a histogram over all the energy 
levels. The HOW method avoids this problem by generating system configurations only in a restricted energy 
interval and estimating the relative weights f2(efc)/fl(ej) of these energy levels from the number of times they 
appear in the sample. From the overlap between energy intervals, one gets the complete fl(efc) function, 
apart from an irrelevant normalization factor. Details about the HOW method can be found in |35|j36|| . 

Finally we note a particularity of this new formalism, which is that the determination of the a constant 
from the normalization condition Pi = 1> require to solve the following equation 

K 

5^n(e fc )F(a + /3e fc ) = l ) (49) 
fc=i 

where efe,/3,f2(efe) are input data for the equation. This equation for a was solved numerically using a 
dicotomic searching method. 

Using the described procedure we have calculated the dependence over the temperature T of the inter- 
nal energy E(N,T), spontaneous magnetization M(N,T) =| ^ i=1 Si |, entropy S(N,T) and free energy 
F(N,T) = E — TS for the 1— dimensional long-range Ising model in the non-extensive regime a = 0.8 and 
the corresponding value x = 1-2 in the Eq. (j|). 

This system has been recently studied within the standard Boltzmann-Gibbs thermostatistics |33| , as well 
as within Tsallis non-extensive ^-formalism In 0] it was numerically verified that the scaling laws of 
the main thermodynamical quantities associated with the Gibbs canonical ensemble are 

E(N,T)/(NN) = e(T/N), 
M(N,T)/N = m(T/N), 
S(N,T)/N = s(T/N), 
F(N,T)/(NN) = f(T/N), (50) 

Since the energy and the entropy scale in different ways, the temperature T has to be scaled as T/N. A 
similar situation arises within Tsallis g-generalized formalism, the concomitant scaling laws being [0] 

E q (N,T)/(NN) = e(TAf(N)/(NN)), 
M q {N,T)/N = m(TAg(N)/ (NN)), 
S q (N,T)/(A q (N)) = S (TA s q (N)/(NNj), 
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F q (N,T)/(NN) = f(TA q (N)/ (AW)), 



(51) 



where 



A q (N) = ( 2 N ^-l)/(l-q), 
A S JN) = (2^1-lVll-gl, 



Af(N) = A q (Nf/A s q (N). 



(52) 



The scaling laws corresponding to the generalized ^-canonical ensemble associated to the weakly nonextensive 
entropy S x are, 

E(N,T)/(NN) = e(TR x (N)), 
M(N,T)/N = m(TR x (N)), 
S X (N,T)/{N^ + N^ X ) = s(TR x (N)), 

F(N,T)/(NN) = f(TR x (N)), (53) 



where 



R X (N) = 



NX + jyl/x 
N*-(2- X )N' 



It is plain from the above equation that 



lim R X {N) = 1, 



so that for large enough values of N the scaling laws ( p3|) become 

E(N,T)/(N*) = e(T), 
M(N,T)/N = m(T), 
S(N,T) X /(N*) = s(T), 
FW,T)/W*) = /(T). 



(54) 



(55) 



(56) 



According to the above equations, the thermodynamic curves of the Ising models ([14|) computed with in- 
creasing N- values must collapse without the need of a temperature rescalation. Numerical evidence of this 
scaling behavior is provided by Figures (||,fj])- It must be realized, however, that the N- values used are not 
large enough to see a complete collapse of the curves depicted. In order to reach complete collapse we need 
R X (N) ps 1 (see the insert of Figure (|b)). 
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FIG. 6. Plot of (a) the internal energy and (b) the magnetization as a function of the scaled temperature. In the 
insert of (b) we plot with solid line the behavior of the scaling factor R X (N) for all the N values and by points for 
the actual simulated sizes N = 100, 400, 1000, 2000, 3000. We use here \ = 1-2. 



We plot in Fig. g the scaled curves for the internal energy and magnetization using the Eqs. (p3|). In the 
insert of that plot we see the asimptotical behavior for the scaling factor R x (N) , and with points we indicate 
the sizes plotted N = 100,400, 1000,2000,3000. Unfortunately we would need very large sizes N ~ 10 s to 
observe an actual scaling without any scale factor, that is R x ~ 1. However from the Fig. 6 one can see that 
there is an actual tendency to the collapse. In fact we expect that the final universal curves will resemble to 
those corresponding to N = 2000, 3000 and we see in the Fig. || a tendency towards a complete collapse. 
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4 8 12 16 



T x R X (N) 

FIG. 7. Plot of the entropy S x as a function of the scaled temperature for the indicated system sizes. In the 
insert we plot the expected no scale temperature behavior. 



Also in Figj7] we show the collapse for the entropy in the main plot with the scaled relations Eqs. ( p3| ) 
and in the insert without scaling. In addition we can see in this plot the actual superadditive behavior of 
the new entropy in an interacting system instead of the simple N independent two states systems of the last 
section. 



6. CONCLUSIONS 

We have shown that the entropy functional S x verifies the main properties usually regarded as essential 
requirements for physically meaningful entropy functionals 0]|,|l2|,[l3| . The entropy S x verifies the first three 
Khinchin axioms. The measure S x satisfies the requirements of positivity, equiprobability, concavity and 
irreversibility. The associated MaxEnt scheme also complies with Jaynes thermodynamic relations. The 
properties verified by S x suggest that it might be a useful measure for describing nonextensive physical 
phenomena as well as for other practical applications. 

The entropic functional S x grows as a power N x of the size N of the system and is thus characterized 
by a weak nonextensive behavior (as opposed to the strong, exponential nonextensivity exhibited by Tsallis 
measure). An attentive reader might think that a power law nonextensive measure can be obtained in a 
trivial way by just defining a new "entropy" equal to (SBoitzmann) x - However, this simple proposal would 
not lead to anything new. The probability distributions maximizing (SBoitzmann) x are the same exponential 
distributions that maximize SBoitzmann- This means that the standard (extensive) thermostatistics would be 
obtained again. On the contrary, the entropic functional S x furnishes a nontrivial new thermostatistics that 
might be appropriate for dealing with some nonextensive systems. As an application of the thermostatistical 
formalism associated with S x we considered a ferromagnetic Ising model with long-range interactions. We 
studied numerically the concomitant scaling laws. 
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Tsallis thermostatistics, which have attracted a great deal of attention recently ||24|— 13 1|] , arise from a 
nonextensive entropic measure S q growing exponentially with the size of the system. Now, if there are 
in Nature extensive systems (that is, those successfully studied within standard statistical mechanics and 
described by Gibbs exponential distributions) as well as Tsallis' strong (exponentially) nonextensive systems, 
it is not unreasonable to expect also the existence of systems exhibiting an intermediate (weak) power-law 
kind of nonextensivity. 

Let us consider a family of dynamical systems characterized by a set of parameters A. We assume that 
there is a certain region S in A-space such that the system's thermodynamic behavior is extensive for A G £ 
and nonextensive otherwise. It would be very surprising if, as a continuous change in the parameters A is 
considered, the system abruptly jumps from an extensive regime to an exponentially nonextensive one. It 
would be physically more sensible if there is an " onset of nonextensivity" boundary between the extensive 
and nonextensive regions in A-space, characterized by a weak, power-law nonextensive regime. A suggestive 
analogy can be made here with the onset of chaos is nonlinear dynamics. Consider, for instance, the well 
known Logistic map. For values of the map's parameter corresponding to chaotic behavior the concomitant 
Lyapunov exponent is positive and near trajectories diverge exponentially in time. However, critical values of 
the map parameter associated with a vanishing Lyapunov exponent still do exhibit a form of "weak chaos". 
In those critical situations, corresponding to the onset chaos, near trajectories do not diverge exponentially 
in time: They diverge in a power-law way p^j3^ |. 

It would be important to gain a further clarification of the physical meaning of the parameter and to 
understand under what circumstances might Nature choose to maximize the functional S x . The only way 
to attain such an understanding is by a detailed study of the dynamics of particular nonextensive systems 
endowed with long range interactions . We hope that our present contribution may stimulate further work 
within this line of inquiry. 
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